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Abstract 

The ^^slow-variation" restriction on the straining flow-velocity gradient field used in the 
Okubo pQ-Weiss [2] criterion is quantified via the Beltrami condition with the divorticity 
framework in 2D hydrodynamic flows. This turns out to provide interesting interpretations 
of the Okubo- Weiss criterion in terms of the topological properties of the underlying vor- 
ticity manifold. These developments are then extended to 2D quasi-geostrophic flows (via 
the potential divorticity framework) and magnetohydrodynamic flows and the Okubo- Weiss 
criteria for these cases are considered. 
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1. Introduction 



A central question in the problem of transport in two-dimensional (2D) turbulent flows 
is how to divide a vorticity fleld into hyperbolic (cascading turbulence) and elliptic (coher- 
ent vortex) regions because the topology of 2D turbulence is parameterized in terms of the 
relative dominance of flow deformation or flow rotation. Okubo [1] and Weiss [2] gave a 
kinematic criterion to serve as a diagnostic tool towards this goal which has been widely 
used in numerical simulations (Brachet et al. |3], Ohkitani Babiano and Provenzale 
[5j) and laboratory experiments (Ouelette and Gollub [6J) of 2D hydrodynamic flows|3 A 
key assumption underlying the Okubo- Weiss criterion is that the vorticity gradient fleld 
evolves quasi-adiabatically with respect to the underlying straining flow-velocity gradient 
fleld. This issue was explored by Basdevant and Philipovitch [9] who tried to improve on 
it by invoking the topological properties of the pressure fleld, while Hua and Klein |10j 
tried to include the strain-rate time evolution explicitly. The purpose of this paper is to 
seek to quantify the ^^slow-variation" restriction on the straining flow-velocity gradient fleld 
used in the Okubo- Weiss criterion via the Beltrami condition with the divorticity framework 
in 2D hydrodynamic flows (Shivamoggi et al. [H]). This also turns out to provide inter- 
esting interpretations of the Okubo- Weiss criterion in terms of the topological properties 
of the underlying vorticity manifold. These developments are then extended to 2D quasi- 
geostrophic flows (via the potential divorticity framework) and magnetohydrodynamic flows 
and the Okubo- Weiss criteria for these cases are considered. 

2. Beltrami Condition Interpretation of the Okubo- Weiss Criterion 

by the following equation 



(la) 
(lb) 



(2a) 
(2b) 



(Ic) 



^The Okubo- Weiss parameter describing the local strain- vorticity balance in the horizontal flow field of a shallow fluid 
layer turns out also to quantify the deviations from two-dimensionality of this flow (Cieslik et al. [7])- More specifically, the 
Okubo- Weiss parameter turns out to be the source turn in the Poisson equation for the pressure (Kamp [8]). 



The vorticity dynamics in 2D hydrodynamic flows is governed 
(Kida [12j, Kuznetsov et al. [13]) 

d'B 

^ = Vx(vxS) 

or 

D'B d'B 

_,_ + (v.V)S = (B.V)v 

where v = {u, v) is the flow velocity, uj is the vorticity, 

a; = V X V 

and B is the divorticity, 

B = V xuj. 

Equation (lb) may be rewritten as 

DB ^ ^ 
-—=A-B 

Dt 
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where A. is the velocity gradient matrix, 



A 



du/dx du/dy 
dvjdx dv/dy 



Sl 
S2 + W 



S2- UJ 



Sl = -2- 



dv 



S2 = 



dv du 
dx dy ' 



u = 



(3) 



dv du 

dy' ^ dx ' dy' dx dy 

If the straining flow velocity gradient tensor Vv is assumed, following Okubo [1] and Weiss 
[2], to temporally evolve slowly so the divorticity field evolves quasi-adiabatically with respect 
to the straining fiow-velocity gradient field, equation (Ic) may be locally approximated by 
an eigenvalue problem with eigenvalues given by 



2 I 2 



UJ 



Q. 



(4) 



The Okubo- Weiss parameter Q is a measure of the relative importance of fiow strain 
(Q > 0, hyperbolic) and vorticity (Q < 0, elliptic). Numerical simulations (Brachet et al. 
[3], Ohkitani Babiano and Provenzale [5]) and laboratory experiments (Ouellette and 
Gollub [6j) of 2D hydrodynamic fiows confirmed that coherent vortices are indeed located in 
elliptic regions while divorticity sheets are locatecH in hyperbolic regions. 

The ^^slow-variation" restriction on the straining flow-velocity gradient field used above 
may be quantified via the Beltrami condition!! with the divorticity framework in 2D hy- 
drodynamic fiows (Shivamoggi et al. [llj). Equations (la-c) yield for the Beltrami state 
(Shivamoggi et al. PJLj), 

'B = av (5) 



a being an arbitrary constant. Using (5), the Okubo- Weiss parameter Q becomes 



/ d^u Y 

\dxdy J 



d'^u d'^u 
dx"^ dy'^ 



(6) 



(6) implies that the Okubo- Weiss parameter also characterizes the topological properties of 
the vorticity manifold - it is in fact the negative of the Gaussian curvature of the vorticity 
surface. Thus, the character of the ensuing time- dependent 2D fiow behavior appears to be 
rooted in the local topological properties of the underlying equilibrium vorticity manifold. 
It may be mentioned that the above reduction was pointed out by Larcheveque [16] on 
the premise of replacing streamlines by isovorticity lines which lacked, as Larcheveque [16] 
admitted, any dynamical meaning - streamlines are actually isomorphic to divorticity lines 
(as implied by the Beltrami condition (5)). 



^It may be noted that divorticity sheets are also more likely to occur near vorticity nulls due to selective rapid viscous decay 
of vorticity in these layers (Shivamoggi et al. |11| '). just as vortex sheets are more likely to form near velocity nulls in 3D 
hydrodynamic flows. 

similar approach was taken previously (Shivamoggi and van Heijst | 14| ) in the quantification of the "slow variation^' 
restriction used in Flierl-Stern- Whitehead jl5) zero angular momentum theorem for localized nonlinear structures in 2D quasi- 
geostrophic flows on the /3-plane. 
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3. The Okubo- Weiss Criterion for Quasi-geostrophic Flows 



Consider a 2D quasi-geostrophic flow in which the barochnic effects are produced by 
the deformed free surface of the ocean. The governing equation (in appropriate units) is 
(Charney [H]) 

§ + (vV)q = (7) 



dt 

where q is the potential vorticity vector 



q = - fc^V + f 



f is the Coriolis parameter, k is the inverse Rossby radius of deformation, k = y fo^ /gH, (/o 
being the local value of |f | and H the depth of the ocean, which is taken to be uniform), and 

V = -V X V (9) 
Upon taking the curl of equation (7), we obtain 

— = V X (v X D) (10a) 

where X> is the potential divorticity vector (in analogy to the potential vorticity vector q), 

D = Vxq = :B + A;V + h (11) 

and 

h = V X f = (/3,0,0) 

/9 being the planetary vorticity gradient. 
Equation (10a) may be rewritten as 

DT) 



Dt 



A ■ D. (10b) 



Following Okubo [Ij and Weiss [2], and assuming that the potential divorticity field evolves 
quasi- adiabatically with respect to the straining fiow velocity gradient field, equation (10b) 
may again be locally approximated by an eigenvalue problem with eigenvalues given by, 

= ^ [UyV^ + Vy^) = Q. (12) 

Equation (10a) yields for the Beltrami state, 

D = 6v (13) 
b being an arbitrary constant. Using (13), the Okubo- Weiss parameter Q becomes 



^-452 



\dxdy J dx"^ dy"^ 



(14) 



which is the same as (6) for 2D hydrodynamic case. This shows that the Okubo-Weiss 
parameter is robust and remains intact under extension to 2D quasi-geostrophic fiows (in the 
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/9-plane approximation to the Coriolis parameter) . The inclusion of the nonhnear variation 
in the Coriohs parameter (the so-called 'y-effect) will, however, lead to changes in the Okubo- 
Weiss parameter. 



dxdy J dx"^ \ dy 



+ 7 



(15) 



4. The Okubo- Weiss Criterion for Magnetohydrodynamic Flows 

Consider a 2D incompressible magnetohydrodynamic (MHD) flow. The equation govern- 
ing the transport of the magnetic fleld B = {Bi, B2) is (Goedbloed and Poedts [I8| ) 

^ = (B • V) V (16a) 

which may be rewritten as 

-— = AB. (16b) 

Dt ^ ' 

If we now assume that the magnetic fleld evolves quasi-adiabatically with respect to the 
straining flow velocity gradient fleld, equation (16b) may again be locally approximated by 
an eigenvalue problem with eigenvalues given by, 

A' = ^ {uyv^ + Vy") = Q. (17) 

Noting that the MHD Beltrami state (Shivamoggi [T^) corresponds to the so-called 
Alfvenic state (Hasegawa [20] ) 

V = cB (18) 
c being an arbitrary constant, (17) becomes 

g = ^ {B,yB2., + 52,/) . (19) 

In terms of the magnetic vector potential A given by 

B = V X A, A = Ai, (20) 

(19) becomes 



Q 



c2 ra^Ay d^Ad^A 

\dxdy J dx"^ dy"^ 



(21) 



(21) implies that the Okubo- Weiss parameter Q for the MHD case characterizes the topo- 
logical properties of the magnetic flux surface - it is the negative of the Gaussian curvature 
of the magnetic flux surface. As with the case of 2D hydrodynamic flows, (21) can serve as 
a useful diagnostic tool to parameterize the magnetic fleld topology in 2D MHD flows. 
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5. Discussion 

The "s/otf variation" restriction on the straining flow-velocity gradient field used in the 
Okubo- Weiss criterion may be quantified via the Beltrami condition with the divorticity 
framework in 2D hydrodynamic fiows (Shivamoggi et al. [H]). This also turns out to 
provide interesting interpretations of the Okubo- Weiss criterion in terms of the topological 
properties of the underlying vorticity manifold. Extension of these considerations to 2D 
quasi-geostrophic fiows (via the potential divorticity framework) shows the robustness of 
the Okubo- Weiss parameter under varying 2D hydrodynamic fiow situations. Extension 
to 2D MHD fiows, on the other hand, provides one again with a useful diagnostic tool to 
parameterize the magnetic field topology in 2D MHD fiows. 
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